Structural anomalies of fluids: 
Origins in second and higher coordination shells 
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Compressing or cooling a fluid typically enhances its static interparticle correlations. However, there are 
notable exceptions. Isothermal compression can reduce the translational order of fluids that exhibit anomalous 
waterlike trends in their thermodynamic and transport properties, while isochoric cooling (or strengthening 
of attractive interactions) can have a similar effect on fluids of particles with short-range attractions. Recent 
simulation studies by Yan et al. [Phys. Rev. E 76, 051201 (2007)] on the former type of system and Krekelberg 
et al. [J. Chem. Phys. 127, 044502 (2007)] on the latter provide examples where such structural anomalies 
can be related to specific changes in second and more distant coordination shells of the radial distribution 
function. Here, we confirm the generality of this microscopic picture through analysis, via molecular simulation 
and integral equation theory, of coordination shell contributions to the two-body excess entropy for several 
related model fluids which incorporate different levels of molecular resolution. The results suggest that integral 
equation theory can be an effective and computationally inexpensive first-pass tool for assessing, based on 
the pair potential alone, whether new model systems are good candidates for exhibiting structural (and hence 
thermodynamic and transport) anomalies. 



I. INTRODUCTION 



A bulk equilibrium fluid is translationally invariant; i.e., its 
one-particle density, p^ 1 \r) = p, is constant. Nonetheless, 
assuming spherically-symmetric interactions, the local den- 
sity pg(r) surrounding a reference particle is a function of 
distance r from its center, where g(r) is the radial distribution 
function (RDF) of the fluid. 1 Although the RDF depends on 
both the form of the interparticle interactions and the thermo- 
dynamic state, some features of its shape are fairly general. 
For example, the RDF vanishes for r less than the effective 
exclusion diameter of the particles. For larger r, it shows an 
oscillatory decay toward unity with peaks loosely correspond- 
ing to coordination "shells". Away from the critical point, the 
structure of the RDF typically persists for distances compa- 
rable to a few particle diameters, reflecting the short range of 
the interparticle correlations. 

Studies of the liquid state have primarily focused on the par- 
ticles in the first coordination shell. This is due in part to the 
important role that nearest neighbors are expected to play in 
determining many physicochemical properties. For example, 
both the non-ideal contribution to the equation of state of the 
hard-sphere fluid (see, e.g., Ref. 2) and the collision frequency 
in Enskog theories for transport processes 3 scale with the 
"contact" density pg{cr), where a is the particle diameter. The 
hard-sphere equation of state is the standard reference system 
for perturbation theories. 1 It also accurately predicts how the 
thermodynamics of "hard-sphere" colloidal suspensions relate 
to their structure, as has recently been experimentally veri- 
fied by confocal microscopy. 4 Furthermore, analysis of first- 
shell contributions to hydration structure and thermodynamics 
helps to understand and make predictions about a wide variety 



of aqueous solution properties ! 5 ' 6 ! 7 ! 8 ' 9 ' 10 ! 11 ! 12 ' 13,14,15 

Although the second shell of the RDF has received com- 
paratively less attention, there is evidence that it also contains 
structural information relevant for understanding nontrivial 
behaviors of liquids. One notable feature is the shoulder— that 
it develops near the freezing transition, which in turn becomes 
a pronounced split pea k 17,18 i 19 ' 20 i 21 i 22 ' 23,24,25 ' 26 in supercooled 
liquid and glassy states. Analysis of the configurations that 
give rise to this structural motif indicate that it reflects frustra- 
tion of icosahedral 23 and emerging crystalline 16 ' 26 order in the 
fluid. Understanding how these these types of structural fea- 
tures connect to relaxation processes of supercooled liquids is 
an active area of research (see, e.g., Refs. l27ll28ll29ll30l) . 

In this work, however, we focus on the second and 
higher coordination shells of the RDF for a different rea- 
son: to understand their role in the structural anomalies 
of fluids. Interparticle correlations of most fluids are en- 
hanced upon (i) compression or (ii) cooling (alternatively, 
strengthening of interparticle attractions). Nonetheless, there 
are a few systems of scientific interest that exhibit no- 
tably different behaviors. For example, compression in- 
duced disordering occurs in water and other fluids with 
anomalous waterlike trends in their thermodynamic and trans- 
port properties, 31 ' 32 ' 33,34 ' 35 ' 36 ' 37 ' 38 ' 39 ' 40 ' 41 ' 42 ' 43 ' 44 Cooling (or at- 
traction) induced disordering, on the other hand, can oc- 
cur in fluids of particles with short-range attractive (SRA) 
interactions, 37 ' 45 such as concentrated suspensions of colloids. 

These anomalies do not appear to be first-shell effects. 
Rather, they reflect how structuring in second and more distant 
coordination shells responds to changes in thermodynamic or 
system parameters. For example, Yan et alM recently demon- 
strated in an insightful paper how the structural anomaly of the 
five-site transferable interaction potential (TIP5P) model^ 6 - for 
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water is quantitatively related to compression induced trans- 
lational disordering of molecules in the second coordination 
shell. Similarly, Krekelberg et al. 45 have shown that the cool- 
ing (or attraction) induced structural anomaly of a square-well 
SRA fluid is due to weakening of second- and higher-shell pair 
correlations. 

The goal here is to study the generality of the above find- 
ings. It is known that a number of models, with vary- 
ing levels of molecular resolution, can qualitatively pre- 
dict the structural anomalies of the aforementioned sys- 
tems, 31 . 32 . 33 - 34 . 35 . 36 ! 37 . 38 . 39 - 40 - 41 . 42 - 43 ! 44 . 45 But do the anomalies 
exhibited by lower resolution models have the same micro- 
scopic origins as those of more detailed models? Moreover, 
can the behavior of the lower resolution models be predicted, 
at least qualitatively, by integral equation theory? If so, it 
would suggest that integral equation theory might serve as a 
valuable first-pass tool in assessing, based on the pair potential 
alone, whether new model systems might be good candidates 
for exhibiting structural anomalies. 

Furthermore, although the structurally anomalous trends 
analyzed here are interesting in their own right, there is a 
more compelling reason to try to understand their origins. 
In short, they appear to be closely linked to other distinc- 
tive dynamic and thermodynamic behaviors. For example, in 
addition to being "structurally anomalous", cold liquid wa- 
ter is also "dynamically anomalous" in that its self-diffusivity 
increases upon isothermal compression and "thermodynam- 
ically anomalous" in that its volume increases upon iso- 
baric cooling. Errington and Debenedetti 31 first noticed that 
these particular anomalies form a cascade in the temperature- 
density plane for the extended simple point charge (SPC/E) 
model 47 of water. 31 Specifically, the thermodynamic anomaly 
occurs only for state points that also exhibit the dynamic 
anomaly. The dynamic anomaly, in turn, is only present for 
states that also exhibit the structural anomaly. Strong correla- 
tions between these three basic types of anomalies have since 
been documented for a wide variety of model systems with 
waterlike properties. 32 - 33 . 36 . 39 - 40 - 43 . 44 - 48 - 49 . 50 - 51 - 52 

A similar connection between structural and dynamic 
anomalies has now also been identified for model SRA 
fluids. 37,45 In those systems, the most commonly studied dy- 
namic anomaly is an increase in self-diffusivity upon cooling 
(or strengthening of interparticle attractions), which can oc- 
cur at sufficiently high particle concentrations! 53 . 54 . 55 . 56 . 57 . 58 . 59 
Krekelberg et a/4 5 - discovered that the self-diffusivity 
anomaly for a square-well SRA fluid occurs only for state 
points that also exhibit the cooling (or attraction) induced 
structural anomaly discussed above. In other words, it ap- 
pears that SRA fluids can also display a cascade of anomalies 
qualitatively similar to those of waterlike fluids. 

Although structural and dynamic properties of these sys- 
tems show unusual dependencies on quantities like temper- 
ature or density, the correlations between structure and dy- 
namics are often similar to those found in simpler liquids 
(e.g., the hard-sphere fluid) . 37 - 38 ' 39 . 45 . 60 In fact, it was recently 
demonstrated 60 that the cascade of anomalies of one water- 
like model system can be semi-quantitatively predicted based 
only on knowledge of the state dependencies of excess en- 



tropy, which measures structural order, 61 and quasi-universal 
excess entropy searings 6 2 . 63 . 64 for the transport coefficients. 
All of this suggests that investigations like the present one, 
which probe the physics of structural anomalies, might also 
provide insights into dynamic and thermodynamic anomalies 
as well. 



II. METHODS 

We used molecular dynamics simulation and integral equa- 
tion theory to examine various models from two classes of flu- 
ids known to exhibit structural anomalies: those with water- 
like properties and those comprising particles with SRA inter- 
actions. For the integral equation theory analysis, we numer- 
ically solved the Ornstein-Zernike equation 65 together with 
an approximate closure relation using the method of Labik 
et al. 66 In the discussion of the models below, we mention the 
specific closures employed and provide further details about 
the molecular simulations. 

We did not perform a systematic study here to determine 
which of many possible closure relations' provides the most 
quantitatively accurate description for each model. Rather, 
our focus was to explore whether integral equation theory 
solved with standard closure relations, such as Percus-Yevick 
(PY) 67 or hypernetted-chain (HNC) 68 , can in fact qualitatively 
predict both the structural anomalies and their microscopic 
origins in the RDF. Molecular simulations of the model sys- 
tems provide the data necessary to make that basic determina- 
tion. 



A. Waterlike fluid models 

We investigated two waterlike models: (1) the SPC/E 47 
model and (2) a lower resolution "core-softened" 36,69 model. 
We chose the SPC/E model because it represents one of the 
most commonly studied effective pair potentials for water, and 
it is known to qualitatively reproduce many of water's dis- 
tinctive thermodynamic, dynamic, and structural properties. 31 
As a result, it provides a reasonable baseline against which 
to compare simpler, lower resolution models. Details of our 
molecular dynamics simulations of the SPC/E model are the 
same as reported in Ref.l37l 

The core-softened mode l 36 - 69 that we studied is more 
schematic. It is defined by the effective pair potential Ucs (?") 
[see Figure [TJa)], 



Ucs(r) = 4e 



-5eexp 



-0.7 



(1) 

where e is the characteristic energy scale. The main idea be- 
hind this potential is that it has two different kinds of repul- 
sions that act at different length scales. The harsh (a/r) 12 
repulsion defines the effective hard-core diameter (<r), while 
the softer Gaussian repulsion extends to considerably larger 
distances. The end result is that the average interparticle 
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separation, and hence the density, of this fluid can depend 
sensitively on both temperature and pressure. The model 
is similar to cold water in that it favors locally open (low- 
density) structures at moderate pressure and low temperature, 
but can collapse to denser structures when compressed or 
heated enough to overcome the soft Gaussian repulsion. Al- 
though this low resolution model does not provide an accurate 
molecular-level description of water, it does qualitatively re- 
produce many of its peculiar thermodynamic, structural, and 
kinetic behaviorsi 36 i 38 - 40 i 69 

To compute the properties of the core-softened model, 
we performed molecular dynamics simulations in the micro- 
canonical ensemble using N = 1000 identical particles of 
mass to. We used the velocity-Verlet technique for inte- 
grating the equations of motion with a time step of At = 
0.002ay/m/e. For the integral equation theory analysis, we 
employed the HNC closure. We chose the HNC approxi- 
mation because of its ability to describe the structure of an- 
other fluid with a soft Gaussian repulsion, the Gaussian-core 
model. 70 We investigated both the SPC/E and core-softened 
models over a wide range of density and temperature, where 
they are known to exhibit structural anomalies ! 31 ' 36 ' 38 ! 40 ' 69 
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FIG. 1: (a) Pair potential of the core-softened model Ucs(r /a)/e 
[see Eq. (Q], (b) Pair potential of the model SRA fluid 
UsRA(r /a) /k&T discussed in the text for various values of polymer 
concentration <f> p . Further details on this SRA model are provided in 
Refs.|Uand|22f 



B. SRA fluid models 

The first SRA fluid model that we considered qualitatively 
describes a solution of (explicit) colloidal particles attracted 
to one another by depletion interactions due to the presence 
of (implicit) non-adsorbing polymers. The details of the col- 
loidal pair potential are provided in Refs. [7l] and[72l but we 



discuss some of its main features below. The colloids are 
spherical and their effective interactions consist of three parts. 
The first is a steeply repulsive, essentially hard-sphere (HS), 
contribution £/hs(?") = ^BT(2a/r) 36 , where 2a is the col- 
loid diameter, k-Q is Boltzmann's constant, and T is tempera- 
ture. The second term represents the aforementioned polymer- 
induced depletion attraction C/aoM' approximated by the 
Asakura-Oosawa 73 potential. The strength of this attraction 
is proportional to the volume fraction of polymers in solution 
(j) p , while the range is controlled by the radius of gyration of 
the polymers R g , set in this case to a/5. A soft repulsion 
C/r is also added to the effective interparticle potential 7 -!- to 
prevent fluid-fluid phase separation. Figure |TJb) displays the 
total colloidal potential Usra = Uns + Uao + Ur for three 
different polymer concentrations <p p . The details of the molec- 
ular dynamics simulations that we performed for this fluid are 
the same as those reported in Ref. U3, with one exception. In 
the original study, a weakly polydisperse system was investi- 
gated. Here, all particles considered had identical radius a and 
mass to. The advantage of focusing on a monodisperse sys- 
tem is that the pair correlations are unambiguously described 
by a single RDF, which facilitates the analysis discussed in 
the next section. For the integral equation theory of this SRA 
fluid, we employed the PY closure. The PY approximation 
is a natural choice here due to its simplicity and its ability to 
describe the structure of liquids with harshly repulsive, short- 
range potentials, 1 (in particular, other SRA fluids 5 - 7 -). 

We also considered a simpler model SRA fluid: a system 
of identical square-well particles with attractive well depth -e 
and width 0.03er, where a represents the hard-core diameter. 
This model is similar to others known to exhibit structural 45 
and dynamic 45 ' 58 anomalies. We also use the PY closure in our 
integral equation theory analysis of this fluid for the reasons 
mentioned above. 



C. Quantification of structural order 

For each of the model fluids, we calculated the state depen- 
dencies of — St/k-B, 



2np / r 2 {g{r)\ng{r) - [g(r) - 1]}, (2) 



£2 



where S2 is the translational pair-correlation contribution^ 5 ^ 
to the excess entropy and p is the number density. We used 
the orientationally averaged oxygen-oxygen RDF in Eq. [2] 
for the analysis of SPC/E water. It has been shown that 
—S2/k-Q not only quantifies the translational order exhibited 
by a fluid (the tendency of pairs of particles to adapt prefer- 
ential separations),— but it also strongly correlates with the 
transport coefficients (see, e.g., Refs. 37,45,64). Other trans- 
lational order parameters have also been introduced to study 
the structure of molecular and colloidal fluids ) 31 ' 61,77 but these 
measures are known to correlate strongly with S2 ) 31 ' 61 and thus 
we exclusively use the latter in our analysis. 

To understand how the various coordination shells of the 
RDF contribute to — S2/&B, we also investigated the cumula- 
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tive order integral I S2 (r), defined as 45 

I S2 (r) = 2irp f r' 2 { 5 (r')W(r') - [g(r') - l]}dr'. (3) 



Note that I S2 (r) — » — ^/^b as r — * oo. 

Finally, we adopted the following criteria to identify struc- 
turally anomalous behavior. 



d[s 2 ] 

Op / rp 

d{- S2 ] 



d[k B T/e] 



< 0, p— structural anomalies (4a) 
> 0, T— structural anomalies (4b) 



As indicated in the Introduction, waterlike fluids ex- 
hibit p— structural anomalies 37 ' 60 and SRA fluids display 
T— structural anomalies^ 



III. STRUCTURAL ANOMALIES 
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A. Waterlike fluids 




1. SPC/E water 

First, we discuss the simulation results for how — S2/fcB 
(i.e., translational order) of the SPC/E water model changes 
with density p. As can be seen in Figure |2ja), SPC/E water 
displays the p— structural anomalies of Eq. (l4ab over the den- 
sity range 0.9g/cm 3 < p < 1.15 g/cm 3 and T < 280K. 
To gain insights into the origins of this behavior, we ex- 
amine the orientationally averaged oxygen-oxygen RDF and 
I S2 as a function of p along the T = 220K isotherm for 
three different density regions: (1) the initial increase of 
—S2{p)/k-Q at low densities [p < 0.9 g/cm 3 , Figs. [2b, c)], (2) 
the anomalous decrease of —S2{p)/k-Q at intermediate densi- 
ties [0.9 g/cm 3 < p < 1.15 g/cm 3 , Figs.HRe)], an d (3) the 
increase of — S2(p)/fcB at high densities [p > 1.15 g/cm 3 , 
Figs.[2f,g)]. 

Compressing the fluid in the lower-density region (1) (p < 
0.9 g/cm 3 ) has relatively little effect on the RDF [Fig. |2b)], 
but it does lead to a small net increase in translational order. 
As can be seen from the behavior of I S2 (r) in Fig. EJc), the 
changes come primarily from the second shell. The reason is 
that the coordination number of water (approximately four, re- 
flecting local tetrahedral hydrogen-bonding to nearest neigh- 
bors) is insensitive to changes in density over this ranged As 
a result, the increase of p is compensated by a slight decrease 
in the first peak of the RDF, and thus the first-shell contribu- 
tion to the structural order remains largely unchanged. In the 
second shell, however, the change in density does not affect 
the RDF (i.e., the strength of the correlations with the central 
molecule). This means that compression induced hydrogen- 
bond bending has allowed more total water molecules into 
the second shell, which in turn leads to an overall increase 
in translational order. 

On the other hand, further increases in density [region (2), 



FIG. 2: Structural data for the SPC/E water model ob- 
tained from molecular dynamics simulations, (a) Structural or- 
der parameter — S2/ks as a function of density p at T = 
220K, 240K, 260K, 280K and 300K. Vertical dotted lines are at 
p = 0.9 g/cm 3 and p = 1.15 g/cm 3 , the approximate bound- 
aries for the region of anomalous structural behavior. (Lower panel) 
Orientationally averaged oxygen-oxygen radial distribution function 
g(r) and cumulative order integral Is 2 ( r ) along the T — 220K 
isotherm [black circles, dashed curve in (a)] for three different den- 
sity regions: (b,c) p < 0.9 g/cm 3 [up to maximum in — S2(p)/&b], 
(d,e) 0.9 g/cm 3 < p < 1.15 g/cm 3 [between maximum and min- 
imum in — S2(p)/&b], (f.g) P > 1.15g/cm 3 [beyond minimum in 
— S2(p)/&b]- The regions are indicated by circled numbers along top 
of (a) and lower panel. In the lower panel, arrows indicate direction 
of increasing density; dashed vertical line is at r = 0.31 nm and dot- 
ted vertical line is at r = 0.57 nm, the approximate locations of the 
first and second minima in g(r), respectively. 



0.9 g/cm 3 < p < 1.15 g/cm 3 ] result in a pronounced de- 
crease in — s 2 (p)/k-Q, i.e. the p-structural anomaly. As can 
be seen in Fig. |2{d), the main implications of compression 
for the interparticle correlations are a dramatic flattening of 
the second coordination shell and an associated shifting in- 
ward of these molecules into the interstitial space between the 
first and second shells. These structural changes are consis- 
tent with the earlier simulation observation 78 that high local 
density can force a fifth molecule from the second shell into 
the periphery of the otherwise four-coordinated first shell. In- 
spection of I S2 (r) [Fig. |2J e )] confirms that the decrease in 
structural ordering is almost entirely due to reduced correla- 
tions between the central and second-shell molecules. In fact, 
Yan et al. convincingly demonstrated that a similar structural 
anomaly in TIP5P water can also be attributed to compression 
induced translational disordering of the second shelli^ 3 - 
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Can these structural changes explain water's self-diffusity 
anomaly? Sciortino et a/. 78 argued, based on molecular simu- 
lation results, that the presence of a fifth molecule in the first 
coordination shell significantly lowers the barriers for trans- 
lational and rotational motions of the central water molecule. 
This suggests that second-shell waters play a central role in 
water's increased mobility under compression. Interestingly, 
since the self-diffusivity of SPC/E water is strongly correlated 
to S2 over these conditions, 37 one can independently draw the 
same conclusion from the data in Fig. [2] 

Finally, we observe that, at sufficiently high densities [re- 
gion (3), p > 1.15 g/cm 3 ], translational order again increases 
upon compression. This is "normal" behavior for dense liq- 
uids, and it simply reflects the fact that smaller volumes 
force particles to adopt locally ordered (i.e., efficient) pack- 
ing structuresi2i&L 



2. Core-softened model 

In this section, we investigate how the translational order of 
the lower resolution, core-softened model of Eq. [T] responds 
to changes in density. First, we consider the results from the 
molecular dynamics simulations. One striking feature of the 
data is that the behavior of — S2/&B as a function of reduced 
density p* = pa 3 , displayed in Fig.Oa), is qualitatively sim- 
ilar to that of SPC/E water [see Fig. |2a)]. Specifically, the 
core-softened model also displays p-structural anomalies over 
the density range 0.08 < p* < 0.175 that become more pro- 
nounced at lower temperature. 

Clearly, the core-softened model is very different from the 
SPC/E model in that the former does not provide a molec- 
ular description of water, and thus it does not favor the for- 
mation of tetrahedrally coordinated hydrogen-bond networks, 
etc. Nonetheless, as we explain below, the main "micro- 
scopic" origins of its density-dependent trends in structural 
order are basically the same as those for the SPC/E model. 

In order to appreciate the similarity between these two 
models, it is helpful to first notice one difference. In the 
SPC/E model, the attractive "hydrogen-bond" interactions 
promote the formation of a first coordination shell, even at 
relatively low density. In contrast, since there are no attrac- 
tions in the core-softened model, the Gaussian repulsion pre- 
vents the "first" coordination shell (near the hard-core diam- 
eter, 1.0 < rja < 1.5) from forming until sufficiently high 
density (p* > 0.1). On the other hand, the "second" coordina- 
tion shell (against the Gaussian repulsion, 1.5 < r/a < 3.5) 
is present even at low density. 

From a qualitative perspective, one might consider each 
core-softened particle as effectively representing a cluster of 
water molecules (e.g., a central water molecule and its four 
nearest neighbors).— In fact, Yan et al. have recently pre- 
sented evidence that a mapping of this sort has quantitative 
merit when one compares, in appropriately reduced form, the 
behaviors of a core-softened ramp model to TIP5P water. 79 
When viewed from this perspective, the formation of the first 
shell in the core-softened model at high density qualitatively 
corresponds, in the molecular picture, to additional water 




FIG. 3: Structural data obtained from molecular dynamics sim- 
ulations of the core-softened potential discussed in the text, (a) 
Structural order parameter — S2/&B as a function of reduced den- 
sity p* = pa 3 at T* = k B T/e = 0.2, 0.3, 0.4, 0.5 and 0.6, where 
a is the particle diameter, and e is the energy scale of the poten- 
tial (see Eq. Q}). Arrow indicates direction of increasing T* , and 
vertical dotted lines are at p* = 0.08 and p* = 0.175, the ap- 
proximate boundaries of the region of anomalous structural behav- 
ior. (Lower panel) Radial distribution function g(r) and cumula- 
tive order integral I S2 (r) along the T* — 0.3 isotherm [red squares, 
dashed curve in (a)] for three density regions: (b,c) p* < 0.08 [up to 
— S2(p*)/kB maximum], (d,e) 0.08 < p* < 0.175 [between maxi- 
mum and minimum in — S2(p*)/A:b], (f,g) p* > 0.175 [beyond min- 
imum in — S2{p*) / ks\- The regions are indicated by circled numbers 
along top of (a) and lower panel. In lower panels, arrows indicate di- 
rection of increasing density; numbers in legends indicate values of 
p*\ vertical dashed line is at r = 1.5er and vertical dotted line is at 
r = 3.5cr, the approximate locations of the first and second minima 
in g(r), respectively. 



molecules (5, 6, etc.) penetrating the first shell of an otherwise 
four-coordinated central water molecule. Once this physical 
relationship between the two models is recognized, the simi- 
larities between their structural properties are easy to under- 
stand. To illustrate this, we carried out a structural analysis of 
the core-softened model identical to that presented above for 
the SPC/E model. 

In particular, we examined the behavior of the RDF and 
I S2 [see Fig. Ob-g)] for the core-softened model as a func- 
tion of density along the T* = k B T/e — 0.3 isotherm 
for three different density regions: (1) the initial increase of 
-s 2 (p*)/k B at low densities [p* < 0.08, Figs.|3b,c)], (2) the 
anomalous decrease of —S2(p*)/kB at intermediate densities 
[0.08 < p* < 0.175, Figs. Hd,e)], and (3) the increase of 
s 2 (p*)/k B at high densities [p* > 0.175, Figs.[3f,g)]. 
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As discussed above, the "first" shell of the RDF is not pop- 
ulated in this model at low density because the core-softened 
particles themselves loosely represent a central water and its 
four nearest neighbors. In this view, the initial compression 
of the core-softened fluid [region (1), p* < 0.08] has an ef- 
fect that is similar to that seen for SPC/E water. The modest 
increase in —S2/k-Q that is observed is due to the increase in 
density and a minor enhancement of structuring in the second 
shell (1.5 < r/a < 3.5). 

Further compression of the core-softened model [region 
(2), 0.08 < p* < 0.175] leads to an anomalous decrease in 
structural order [Fig. a)]. Figs. |3jd) and (e) indicate that 
the disordering is again due to a flattening and shifting inward 
of the second shell. Moreover, the "first" shell of the core- 
softened particles begins to emerge, which schematically rep- 
resents, in the approximate molecular view discussed above, 
that additional water molecules are effectively penetrating into 
the four-coordinated first shell. 

Similar to SPC/E water, it is known that there is a strong 
correlation between excess entropy and self-diffusivity for the 
core-softened model. 38 This information, together with the re- 
sults shown here, support the view that the self-diffusivity 
anomaly of the core-softened model is also linked to its 
density-dependent second-shell structure. 

As expected, at higher density [region (3), p* > 0.175], 
compression leads to an increase in structural order due to 
simple-liquid-like structuring of particles in the first coordi- 
nation shell [Figs. [3jf,g)] . In short, the qualitative response to 
changes in density of the structural order and its coordination- 
shell contributions for the core-softened model are very sim- 
ilar to those of the more detailed SPC/E water model. This 
finding is consistent with the recent demonstration that one 
can approximately map the anomalies of TIP5P water onto 
those of a similar two-scale ramp potential^ 

As a final point, we show in Fig.|4]that the integral equation 
theory of the core-softened model can qualitatively predict all 
of the trends shown in Fig. [3] The ability of this approach 
to reproduce the structural features seen in simulations, to- 
gether with the quasi-universal connection between structure 
and transport coefficients of liquids ^ 62,63 ' 64 suggests that inte- 
gral equation theory might serve as a valuable first-pass tool in 
assessing whether other model systems represent good candi- 
dates for exhibiting static and dynamic anomalies. However, 
if the intention is to ultimately use it as a quantitatively accu- 
rate predictive tool, then more comprehensive investigations 
of alternative closure relationships, in the spirit of Ref. |3Jj 
will be necessary. 



B. SRA fluids 

One of the key aspects of short-range attractive (SRA) flu- 
ids is that their structurally anomalous behavior occurs as a 
function of the reduced interparticle attractive strength e/k^T 
at constant particle packing fraction <fr c , where — e represents 
the well depth of the interparticle attraction. In most typi- 
cal atomic or molecular fluids, one finds that structural order 
(— S2/fce) increases with tjk^T. SRA fluids are anomalous in 




FIG. 4: Structural data for the core-softened waterlike model from 
integral equation theory, (a) Structural order parameter —silks, 
as a function of reduced density p* — pa 3 at the same values of 
T* = fer/t as in Fig. [3J a )> where a is the particle diameter, 
and e is the energy scale of the potential (see Eq. QJ). Arrow in- 
dicates direction of increasing T* , and vertical dotted lines are at 
p* = 0.075 and p* = 0.165, the approximate boundaries of the 
region of anomalous structural behavior. (Lower panel) Radial dis- 
tribution function g(r) and cumulative order integral I 32 (r) along 
the T* = 0.3 isotherm [red dashed curve in (a)] for three different 
density regions: (b,c) p* < 0.075 [up to — S2(p*)//cb maximum], 
(d,e) 0.075 < p* < 0.165 [between maximum and minimum in 
— S2(p*)/A:b], (f,g) p* > 0.165 [beyond minimum in — S2(p*)/te]. 
The regions are indicated by circled numbers along top of (a) and 
lower panel. In lower panels, arrows indicate direction of increasing 
density; numbers in legends indicate values of p*\ vertical dashed 
line is at r — 1.5a and vertical dotted line is at r = 3.5<r, the ap- 
proximate locations of the first and second minima in g(r), respec- 
tively. 



that, at sufficiently high values of <j> c , the opposite trend can 
be observedri 4 ^ i.e., attractions counterintuitively decrease 
the amount of structural order. 

In this section, we briefly discuss how we used molecular 
simulation and integral equation theory to gain insights into 
this trend. We accomplished this by exploring the various 
coordination-shell contributions to — S2/ &b for the two model 
SRA fluids discussed in Section llTBl 



1. Colloid-polymer mixture 

We begin by investigating the behavior of the model 
colloid-polymer systemic by molecular simulation. The ef- 
fective colloid-colloid pair potential for this model was pre- 
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FIG. 5: Structural data obtained from molecular dynamics simula- 
tions of the model colloid-polymer SRA fluid discussed in the text, 
(a) Structural order parameter — sz/ks as a function of polymer vol- 
ume fraction cj> p (i.e., strength of colloid attractions) at colloid pack- 
ing fraction cj> c = 0.4. Vertical dotted line at <j> v = 0.1, the location 
of the minimum in —S2(<j> p )/kB- (Lower panel) Radial distribution 
function g(r) and cumulative order integral I S2 (r) along the isochore 
4> c — 0.4 [black circles in (a)] for two polymer concentration ranges: 
(b,c) (f>p < 0.1 [below minimum in — S2(cA p )/fcB], (d,e) cA p > 0.1 
[above — S2{<f> P )/k-B minimum]. The regions are indicated by circled 
numbers along top of (a) and lower panel. In lower panels, arrows in- 
dicate direction of increasing cj> p \ the parameter a indicates colloidal 
particle radius; vertical dashed line is at r — 3a and vertical dotted 
line is at r — 5a, the approximate locations of the first and second 
minima in g(r), respectively. 



sented earlier in Fig. 02b) for several values of polymer pack- 
ing fraction <fi p . Since the reduced well-depth of this potential, 
e/k-oT, scales as <f> p , we analyze structural order below as a 
function of the latter. 

In particular, Fig. [5J a ) illustrates how — S2/fce varies as a 
function of cj) p at a particle packing fraction of cf) c — 0.4. As 
expected for SRA fluids, — S2{4> p )/k-Q exhibits a minimum at 
(j> p « 0.1. In other words, this fluid displays the structural 
anomaly of Eq. d4bl for cf> p < 0.1. To understand the origins 
of this trend, we study the RDF and the cumulative order in- 
tegral I S2 (t) as a function of <p p in two qualitatively different 
regions: (1) the anomalous decrease in — S2((j) p ) / k-Q at low 
polymer concentrations (low interparticle attractions), and (2) 
the "normal" increase of — S2 (4> p ) /fee at higher polymer con- 
centrations (higher interparticle attractions). 

What specific changes to coordination shell structure ex- 
plain the attraction induced disordering that occurs at small 
4> p [region (1), <fi p < 0.1]? First, note that strengthening in- 



terparticle attractions considerably increases but narrows the 
first peak of the RDF [Fig.|5fb)]. These two effects essentially 
cancel so that the first-shell contributions to I S2 (r) are insensi- 
tive to <f> p over this range [see Fig. [3c)]. However, attractions 
also slightly shift the higher coordination shells of the RDF 
inward and diminish their overall correlation with the central 
particle. These latter modifications to the structure of the sec- 
ond and higher coordination shells give rise to the anomalous 
decrease in the structural order of this system. They are also 
consistent with behavior observed in the recent Krekelberg et 
a/4 5 - simulations of the square-well SRA fluid discussed in the 
Introduction. A microscopic interpretation of this trend is that 
SRA interactions drive weak particle clustering at low <f) p (ex- 
plaining the sharpening and narrowing of the first peak). This 
clustering, in turn, opens up channels of free volume in the 
fluid and disrupts the uniform hard-sphere-like packing order 
in the second and higher coordination shells . 59 ' 74 ' 80 ' 81 

Under conditions where the aforementioned structural 
anomaly occurs, increases in <j) p also increase the mobility 
of the fluid. 45 Very similar to the waterlike fluids discussed 
above, it is known that S2 and self-diffusivity are strongly cor- 
related for the model colloid-particle mixture^ 7 - As a result, 
the self-diffusvity anomaly appears to also derive from sub- 
tle structuring effects in the second and higher coordination 
shells. 

As one would expect, however, increasing <f> p ultimately in- 
creases structural order, if the interactions are sufficiently at- 
tractive [region (2), <\> p > 0.1]. The attractions lead to the 



formation of strongly bonded particle clusters 



59.74.80.81 



which 



is reflected by the increased height of the first peak of the RDF 
[Fig-Ed)] and the associated rise in the first-shell contribution 
to/ S2 (r) [Fig. He)]. 

In closing, we test in Figs. [6] and [7] whether integral equa- 
tion theory is able to qualitatively capture these attraction in- 
duced structural changes for both model systems introduced 
in Section |HB] the colloid-polymer fluid and the square-well 
fluid, respectively. Comparison of Fig.|6]with Fig.|5]and Fig. [7] 
with Figure 4 of Ref. 45 demonstrate that this is indeed the 
case. The success of these predictions strengthens the case 
that integral equation theory will be a useful tool in assessing 
whether future model systems of interest might display struc- 
tural anomalies. 



C. Conclusions 

Although the structural order of a fluid is usually enhanced 
by isothermal compression or isochoric cooling, a few notable 
systems show the opposite behaviors. Specifically, increas- 
ing density can disrupt the structure of waterlike fluids, while 
lowering temperature (or strengthening of attractive interac- 
tions) can weaken the correlations of fluids with short-range 
attractions. The two-body translational contribution to the ex- 
cess entropy provides a quantitative measure of these changes 
in structural order. It is a particularly insightful quantity to 
study because (i) its contributions from the various coordina- 
tion shells of the radial distribution function can be readily 
determined, and (ii) it correlates strongly with self-diffusivity, 
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FIG. 6: Structural data for the model colloid-polymer SRA 
fluid discussed in the text from integral equation theory. (a) 
Structural order parameter —S2/kB as a function of poly- 
mer volume fraction <j} p at colloid packing fractions (j> c = 
0.3, 0.325, 0.35, 0.375, 0.4, 0.425, 0.45, 0.475 and 0.5. Arrow in- 
dicates direction of increasing cj> p , and vertical dotted line is at 
4>p = 0.1, the approximate boundary of the region of anomalous 
structural behavior. (Lower panel) Radial distribution function g(r) 
and cumulative order integral I S2 (r) along the isochore <j> c = 0.475 
[dashed violet curve in (a)] for two polymer concentration ranges: 
(b,c) <f)p < 0.1 [below minimum in — S2(<f} p ) / ks], (d,e) (f> p > 0.1 
[above ~S2(4> P ) /ks minimum]. The regions are indicated by circled 
numbers along top of (a) and lower panel. The parameter a indicates 
colloid radius. In lower panels, arrows indicate direction of increas- 



ing ( 



, vertical dashed line is at r — 3a, and vertical dotted line is at 



r = 5a, the approximate locations of the first and second minima in 
g(r), respectively. 
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FIG. 7: Structural data for the square-well fluid discussed in the text 
obtained from integral equation theory, (a) Structural order param- 
eter — S2/fcB as a function of reduced attractive strength e/ksT at 
particle packing fractions <f> c — 0.4, 0.45, 0.5, 0.525, 0.55, 0.56 and 
0.57. Arrow indicates direction of increasing <j) c , and vertical dot- 
ted line is at e/fceT = 0.9, the approximate boundary of the region 
of anomalous structural behavior. (Lower panel) Radial distribution 
function g(r) and cumulative order integral I 32 (r) along the iso- 
chore 4> c = 0.55 [dashed blue curve in (a)] for two attractive strength 
ranges: (b,c) e/ksT < 0.9 [below minimum in —S2(fi/kBT)/kn\, 
(d,e) e/ksT > 0.9 [above — S2(e/fcBT)/fcB minimum]. The re- 
gions are indicated by circled numbers along top of (a) and lower 
panel. The parameter a indicates colloid diameter. In lower panels, 
arrows indicate direction of increasing e/ksT; vertical dashed line is 
at r = 1.4cr and vertical dotted line is at r = 2.3<r, the approximate 
locations of the first and second minima in g(r), respectively. 



which allows it to provide insights into the dynamic anomalies 
of these fluids. 

Here, we have presented a comprehensive study, by both 
molecular simulation and integral equation theory, of the co- 
ordination shell contributions to the two-body excess entropy 
for several model systems. These models incorporate different 
levels of molecular resolution, but all exhibit the aforemen- 
tioned structural anomalies. The results of this study support 
the emerging view that the structural anomalies of these fluids 
can generally be attributed to quantifiable changes in the sec- 
ond and higher coordination shells of the radial distribution 
function. They also demonstrate that integral equation theory 
can serve as an effective first-pass tool for assessing, based on 
the pair potential alone, whether new model systems are good 
candidates for exhibiting static and dynamic anomalies. 
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